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1. Introduction 

Investigations of the phase structure of strongly-interacting matter have received a consider- 
able amount of attention in the last years, both, theoretically as well as experimentally. Among the 
most prominent features of the QCD phase diagram is the crossover from a confined phase with 
spontaneously broken chiral symmetry to a chirally symmetric and deconfined phase. A confined 
phase can be linked to a ground state that respects center symmetry, 1 and the associated order pa- 
rameter is the Polyakov loop. It vanishes in the center symmetric phase and becomes finite as soon 
as the center symmetry is broken [1]. At vanishing quark chemical potential both transitions occur 
roughly at the same temperature which led to the idea and introduction of new dual observables in 
lattice QCD [3-6]. In particular, the dual chiral condensate pioneered in [3] is constructed from 
the chiral condensate, the order parameter of chiral symmetry breaking. More recently, dual ob- 
servable have also become accessible within functional methods [7-9] and have been successfully 
applied to investigate the center transition. In the present work novel order parameters for the center 
symmetry and its breaking are introduced and analyzed in QCD as well as fundamentally charged 
scalar QCD. The order parameters are determined by the corresponding matter propagators without 
any additional renormalization. 

2. (Scalar) Quantum Chromodynamics 

We address the deconnnement transition in ordinary QCD as well as in scalar QCD, where the 
quarks are replaced by fundamentally charged scalars (see e.g. [10]) both in Landau gauge. 

The matter propagators are calculated by means of the corresponding Dyson-Schwinger equa- 
tion (DSE) [11]. As an example, the DSE for the quark propagator is shown diagrammatically in 
Fig. 1 , where thin lines and dots represent bare propagators and one-particle irreducible vertices 
while thick lines and dots denote the corresponding dressed quantities. Accordingly, the DSE for 
the scalar propagator, which is shown in Fig. 2, is more involved and contains more diagrams due 
to the presence of additional bare vertices such as the scalar self-interaction and the quartic scalar- 
gluon vertices. In one-loop approximation only the momentum-independent tadpole diagrams are 
left in addition to the gluon exchange diagram. The tadpoles, however, can be treated by adjusting 
the renormalization constants appropriately. Hence, in a one-loop approximation, the DSE for the 
scalar propagator is of the same structure as the one for the quark propagator, cf. Fig. 1. 

Explicitly, at finite temperature T the DSE for the quark propagator S(p) reads 

S-\p)=Z 2 S,\p)-Z lF C F g 2 T £ f^- 3 f S (k)T v (k,p; q )D^(q) (2.1) 



and correspondingly for the scalar propagator D${p) 




In a strict sense center symmetry is realized only in the limit of infinitely heavy quarks while in real QCD the 
symmetry is always explicitly broken, see e.g. [1,2]. 



2 



Center Phase Transition from Matter Propagators 



Mario Mitter 



Figure 1: DSE for the quark propagator. 



-i 



2 




i 



2 




Figure 2: DSE for the fundamentally charged scalar propagator. 



For the four-momenta we use k = (k,(Oi c (d)) and the gluon momentum is constrained by momentum 



conservation to q = p — k. The wave function renormalization of the quark (scalar) fields are 
denoted by Z2 (Z2), the renormalization constants of the quark-gluon (scalar-gluon) vertex by Z\p 
(Z\p) and Z m labels the scalar mass renormalization constant. The quadratic Casimir invariant in 
the fundamental representation of the gauge group SU (3) is Cf = 4/3 and the coupling constant at 
the renormalization scale is given by g. 

In general, exp(/0)-valued boundary conditions in the fourth spacetime direction are realized 
by introducing generalized Matsubara frequencies (0^(6) = (27ini ( + d)T, where the sums over 
(Qk(0) in the DSEs run over the corresponding discrete values G Z. The usual (anti-)periodic 
boundary conditions for (fermions) bosons are obtained by setting (6 = n) 8 = 0, respectively. 

Both DSEs depend on the gluon propagator D^ v as well as on the corresponding matter-gluon 
vertices T v and Fj. Solutions of the DSE for the Landau gauge gluon propagator at non- vanishing 
temperatures have been obtained in [12]. In addition, data are available from (quenched) lattice 
simulations at finite temperatures and have already been successfully implemented in functional 
equations for the quark propagator [8, 13] 2 . In this work we will apply the fit functions for the 
gluon propagator proposed in [13] and hence we omit explicit expressions. 

For the matter-gluon vertices the situation is less satisfactory since the temperature behavior 
of these vertices is not known generally. Some modeling has to be employed such as in [8] where 
the following expression 



for the quark-gluon vertex can be found and will be used also in this work. The Ansatz is motivated 
by Slavnov-Taylor identities of the Abelian gauge theory (see e.g. [15]) and by the running of the 
non-perturbative coupling of the Yang-Mills theory. The purely phenomenological parameters di 

2 Recently, also unquenched lattice data for the Landau gauge gluon propagator are available [14], which will be 
used in future studies of the system. 




(2.3) 
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Figure 3: Solution of the truncated DSE Eq. (2.2) for the scalar propagator as function of p 2 = x at the 
lowest Matsubara frequency for periodic 0=0 (left panel) and antiperiodic 9 = % (right panel) boundary 
conditions at different temperatures with renormalization scale ji = 4 Ge V and mass m =1.5 GeV. 



and are specified in [13], whereas the anomalous dimension 28 = —18/44 and j8o = l\N c /3 
ensure a correct perturbative running coupling in the UV regime for SU (N c ) gauge theory. The 
renormalization scale of the Yang-Mills sector has been fixed by a(ju) = 0.3 and A = 1.4 GeV. 
The factor Z3 allows to apply the Slavnov-Taylor identity Z\p = Z2/Z3 in Landau gauge [16] which 
yields finally only a dependence on the quark wave function renormalization Z2. The gluon mo- 
mentum is denoted by q and k and p are the in- and outgoing quark momenta, respectively. 

In a similar context we employ 

x - D^ l (p 2 )-D^(k 2 ) , 
r v s (k,p;q)=Z 3 s {I p i_ k2 S { ] (p + k) v 

A A 2 q 2 ^oa(M)ln[?7A 2 + l] V g l ^ 

XdY \A 2 + q 2+ q 2 +A 2 { 4x ) } 

for the scalar-gluon vertex, where one additional parameter d\ = 0.53 has been introduced and all 
remaining parameters are the same as in the quark-gluon vertex. In contrast to the quark-gluon 
vertex, the scalar-gluon vertex is dressed only with the vacuum propagators D^ l (p 2 ). For the nu- 
merical solution of the corresponding DSEs we rewrite the propagators S^ 1 (p) = iy^(O p (d)C{p) + 
ipA{p)+B(p) and D${p) =Zs(p 2 ,C0 p (8))/(p 2 + (O p (d) 2 ) in terms of dressing functions in a stan- 
dard way. Details on the numerical implementation as well as on the renormalization scheme will 
be published elsewhere [17], cf. also [18, 19]. 

Numerical results for the scalar propagator around temperatures of the center transition in 
the quenched theory with T c « 277 MeV are shown in Fig. 3 for periodic (left panel) as well as 
antiperiodic boundary conditions (right panel). The mass of the scalars has been fixed to m = 1.5 
GeV which results in a quite inert behavior around T c . These results demonstrate that there are no 
direct modifications of the scalar propagator in the vicinity of the transition. This motivates the 
construction of more sensitive order parameters for the center phase transition. 
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3. Center Phase Transition and Dual Order Parameters 



Order parameters for the center phase transition can be composed with functional methods 
via dual observables like e.g. the dual chiral condensate [7, 8] or the dual density [9]. As has 
been discussed in the latter reference, dual quantities can be evaluated in two different ways. In 
general, dual order parameters are constructed from some boundary-condition dependent operator 
Oq, where 6 denotes the phase of the U (l)-valued boundary conditions. Originally, such operators 
were introduced in lattice calculations [3-6] and evaluated in QCD with the original boundary 
conditions, i.e., with (anti-)periodic boundary conditions for (fermions) bosons. 

On the other hand, these operators can also be evaluated with functional methods in various 
theories with general boundary conditions, referred to as QCDg in [9]. However, dual observables 
evaluated in this way can serve as order parameters only if the deconfinement transition temperature 
at physical boundary conditions is a lower bound for the transition temperatures in the different 
theories QCD e [9]. 

In previous studies of the center phase transition the dual chiral condensate has been calculated 
with functional methods based on the quark propagator [7, 8, 13]. It is given by an expansion in 
complex Fourier modes 

fin aq 

Z„ = J o ^e- M {W)B (3-D 



with a 6 -dependent quark condensate 



(W)e=Z 2 N c T £ f ^P t r D S(p,(O p (e)) . (3.2) 



a p {ey 



(2k) 



For arbitrary n that is not a multiple of N c , £„ can then serve as an order parameter for center 
symmetry. Usually, the dual chiral condensate Ei is used, also called dressed Polyakov loop. 
In a lattice formulation it contains contributions from all time-like loops around the torus with 
winding number n = 1 [3, 4, 6] and transforms similar to the ordinary Polyakov loop [1] under 
center transformations. 

Hence, the calculation of the dual chiral condensate requires the chiral condensate with general 
boundary conditions, which has to be regularized for non-vanishing quark masses. As an already 
finite alternative we propose 



old " 1 : 



■'0 z/t fa- 



l -tr D S(Q,(O p {e)) 



(3-3) 



a p (6) 

as an order parameter for the center phase transition. Zq q is finite, because the sum scales like 
I /(Op for large Matsubara modes. Similarly, for scalar QCD we propose the order parameter 

1-271 JQ 

= / ^-e-' e Z Sie , Z 5 , e = T £ D 2 s (0,co p (d)) . (3.4) 
More details on these order parameters will be presented in an upcoming publication [17]. 
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Figure 4: Left: g of Eq. (3.4) as a function of the boundary conditions for temperatures around the 
transition. Right: The dual condensate E5 as a function of the temperature {pi = 4 GeV and m = 1.5 GeV). 
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Figure 5: Parameter dependency of the dual condensate E5 as a function of temperature for three different 
values of the vertex parameter d\ (left) and for three different values of the mass of the scalar field (right); 
(ju = 4 GeV and m = 1.5 GeV). 



4. Results 

In order to confirm that both quantities, Eg in Eq. (3.3) and E5 in Eq. (3.4), are well-defined 
order parameters we investigate both theories, QCD and scalar QCD, at finite temperatures. 

In the left panel of Fig. 4 the -dependence of JLg,e f° r the quenched scalar theory is shown for 
temperatures around the transition. From the definition Eq. (3.4) it is clear that £5 vanishes as long 
as the L s q is constant while a 6 -dependency is necessary for a non- vanishing order parameter. 
However, these findings for L s q are qualitatively similar to QCD with finite quark masses [3]. 
In the right panel of Fig. 4 the dual condensate £5 is shown as a function of the temperature 
which nicely demonstrates its property as an order parameter for the center symmetry. Below the 
transition temperature of the quenched theory around T c « 277 MeV it vanishes and is finite at 
higher temperatures. In the vicinity of the critical temperature the temperature behavior of the 
order parameter depends crucially on the model parameters used for the scalar-gluon vertex. This 
is demonstrated in Fig. 5, where the temperature behavior of £5 is shown for different values of the 
d\ parameter in the scalar-gluon vertex Eq. (2.4) (left panel) and for various values of the mass of 
the scalar field (right panel). Stronger deviations from a vanishing order parameter slightly below 
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8/271 T [GeV] 



Figure 6: Left panel: The quark condensates (y/y/)e and Eg g, as defined in Eq. (3.3), as a function of 
the boundary angle for different temperatures in the chiral limit (dashed lines: T = 273 MeV, solid lines: 
T = 283 MeV). Right panel: The order parameters Eg and Ei as defined in Eq. (3.1) and (3.3) as a function 
of the temperature. 

the transition temperature are observed for smaller mass. This indicates that the vertex sensitivity 
increases towards smaller masses. 

Finally, the proposed new order parameter for ordinary QCD shows a similar behavior which 
is presented in the left panel of Fig. 6. In the figure the 0-dependency of Lqq in comparison to the 
chiral condensate is plotted as a function of the generalized boundary conditions for two different 
temperatures. The results have been obtained in the chiral limit of the quenched theory. Above 
the transition temperature, Lq q shows the same characteristic plateau as the chiral condensate for 
0's close to the physical antiperiodic boundary conditions 6 = %. Due to the restoration of chiral 
symmetry the condensate vanishes which also affects Lqq through the scalar dressing function 
B(p) in the quark propagator. The slight variations of Lq in the center symmetric phase below T c 
results from a non-constant Lq q and can be attributed to lattice artifacts as well as the choice of 
the quark-gluon vertex model. However, this effect is more pronounced in Lq^q than in the chiral 
condensate. The right panel of Fig. 6 shows the corresponding order parameter Lq in comparison 
to the dual chiral condensate. Both vanish below T c and jump immediately to non- vanishing values 
above T c . Their deviation at larger temperatures can be assigned to different dimensionalities. 



5. Conclusions 

We investigated the center phase transition of QCD as well as fundamentally charged scalar 
QCD in a quenched formulation. Novel order parameters are proposed along the lines of previously 
constructed dual observables accessible by functional methods. Solving the Dyson-Schwinger 
equations for the corresponding matter propagators numerical results for these order parameters 
are presented. A parameter dependency of the employed matter-gluon vertices on the presented 
results is found and motivates a more detailed investigation of its temperature dependence (see 
also [20] for corresponding investigations at vanishing temperature). 
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